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ELECTROMAGNETIC WAVE PROPAGATION IN A COLD, 

COLLISIONLESS ATOMIC HYDROGEN PLASMA 
by Richard R. Woollett 
Levis Research Center 


SUMMARY 

The characteristics of electromagnetic vave propagation in a cold, colli- 
sionless plasma immersed in a steady, uniform magnetic field are presented in 
dimensional units. Data for the "principal” waves are presented as curves of 
vave number, vavelength, phase velocity, and group velocity as functions of fre- 
quency. Wave numbers for oblique propagation are presented for angles of 0.1°, 
1°, 10°, 45°, 80°, 89°, and 89.9° relative to the magnetic-field direction. The 
range of magnetic-field strengths covered by the plots is 10 2 to 10 5 gauss, and 
the range of electron densities is 10 10 to 10 15 ions per cubic centimeter. Cut- 
off and resonance frequencies are also plotted as functions of density and 
magnetic-field strength. 


INTRODUCTION 

The intense interest shovn In plasma vaves has led to several independent, 
though similar, theoretical treatments of the problem. Results have appeared in 
a variety of forms. The case of vave propagation in a cold, collisionless 
plasma is typical in this regard. Although the analysis of the cold plasma has 
been quite exhaustive, the vave characteristics have been presented either as 
sketches shoving the limiting cases in detail or as detailed plots covering re- 
stricted ranges of the propagation frequency. These presentations used nondi- 
mensional expressions over a restricted frequency region in order to correlate 
the data from many conditions on a single plot. The experimentalist, hovever, 
usually prefers dimensional data that cover a broad range of application. In 
this report the more common vave propagation parameters are presented as dimen- 
sional quantities and are plotted over an extended frequency range. These 
curves are intended to complement the nondimen sional plots that have been pre- 
sented over various frequency ranges (refs. 1 and 2). 

Characteristics of the plasma vave commonly used in microvave diagnostics 
to measure electron density (termed the ordinary vave) are independent of mag- 
netic field. Therefore, its vave properties can be presented conveniently in 
either nondimensional or dimensional form, because the only parameter is electron 
density. It is possible, hovever, to devise diagnostic techniques based on char- 
acteristics of other types of vaves. Since the plasma densities that can be de- 
termined from these other vaves generally depend on magnetic field, care is re- 
quired In selecting the propagation frequency. This detailed selection of test 
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conditions is made more convenient "by the use of comprehensive plots that cover 
a large range of experimental frequencies. 

The curves presented in this report apply for small- amplitude electromag- 
netic waves (since a perturbation technique is employed) propagating through an 
atomic hydrogen plasma immersed in a uniform steady magnetic field. Calculations 
made from Astrom 1 s equations (ref. 2) predict, for a cold, collisionless plasma, 
all the important electromagnetic wave modes and their resonances. Included in 
the results are curves of wave number, wavelength, phase velocity, and group ve- 
locity of waves propagating parallel and normal to the magnetic field over a 
spectrum of excitation from low-frequency radio waves (1.6x10 s cps) to submilli- 
meter microwaves (l. 6xl0 12 cps). In addition, wave-number plots for propagation 
at angles other than 0 o and 90° are presented. Ranges of magnetic-field inten- 
sity and particle density covered by the graphs are, respectively, 10 2 to 10 5 
gauss and 10 10 to 10 15 ions per cubic centimeter. 


THEORY 

The basic equations for the propagation of an electromagnetic wave in a 
cold plasma with overall charge neutrality have been derived from various mathe- 
matical models. The contribution that is usually accepted as the original major 
work is that by Astrom (ref. 2). He assumed an infinite, homogeneous medium of 
ions and electrons that were motionless in the unperturbed state and were im- 
mersed in a uniform static magnetic field. The plasma was considered collision- 
free. Therefore, the only forces to act on the charged particles were those 
originating from the macroscopic electric and magnetic fields. Since there were 
assumed to be no collisions between particles, the presence of neutral constitu- 
ents would not alter the results obtained from this theory. The ion and electron 
velocities were obtained by solving the equation of motion for a charged particle 
acted on by the Lorentz force. The plasma currents obtained from these particle 
velocities were substituted into Maxwell’s equations to obtain a self-consistent 
solution in terms of the dielectric constants of the medium. Although Astrom as- 
sumed a plane wave in his development of the dielectric coefficients, Epstein 
(ref. 3) obtained the same dispersion relation for cylindrical waves. 


Dispersion Relation 

There are various forms in which the resulting equations may be presented, 
and each form is more convenient under certain circumstances than the others. 

The equations considered in this report involve the various components of a modi- 
fied dielectric tensor K; the modified dielectric coefficients used herein are 
obtained by multiplying Astrom 1 s dielectric coefficients by (cjd/c) 2 . Components 
of the tensor are K^, K^, Kp, and K s , where the subscripts 7, r, p, and s 
refer to left rotation, right rotation, axial, and transverse electric-field 
polarizations. (Symbols are defined in appendix A. ) The equations used in this 
report can be derived from the dispersion relation (ref. 4) 

j/k 4 - #k 2 +< 4=0 (1) 
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where 


st = Kg *sin 2 0 + Kp cos 2 0 
& = 1^2 sin 2 0 + KpK g (l + cos 2 0) 
<€ -K pKr K 2 




It should he noted that Kp is independent of magnet ic-f ield strength. 

Waves propagating at 0° and 90° with respect to the magnetic field are re- 
ferred to as principal waves or modes. Waves propagating at angles between 0° 
and 90° are referred to as oblique waves. If a rectangular coordinate system is_ 
selected so that the positive directions of x and _z are parallel to (Bxk) xB 
and B, respectively, the coordinate components of k become 

y ( 5 ) 

where the numbers in the subscript identify the two roots obtained from equa- 
tion (l). The components k x and k z , which are most commonly used in bound 
plasma problems, are the quantities that are generally measured experimentally. 

Another convenient form in which the propagation vector may be presented is 
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( 6 ) 


tan 2 © = 


- S ( k2 ~ K r )( k2 - K z ) 

(k 2 - Kp)(K s k 2 - Kr K 2 ) 


This form immediately yields the wave numbers of the principal plasma modes by 
setting 0=0° and 90°. Hence, 



for k x = 0, 6 = 0 C 


for k z = 0, 6 = 90 c 


(7) 


where k 2 and k^. are the wave numbers of the left and right circularly polar- 
ized plane waves propagating in the direction of the magnetic field, kp is the 

wave number of the P mode (ordinary wave) in which the electric vector is par- 
allel to the field direction, and k s is the wave number of the S mode (extra- 
ordinary wave) in which the electric vector is transverse to the magnetic-field 
direction. The only dispersion relation of the principal modes that is not di- 
rectly represented in equation (7) is the one for the S mode, which can be 
written separately as 

k 2 _ o£ + ^ci 03 ~ ^ce^ 0 ~ ^ci^ce ~ ^p)^ 2 ~ ^ci^ + ^ce^ - ^ci^ce “ ^p) 
c 2 (cd 4 - O^CD 2 - CD 2 e 60 2 + co 2 ^ - o|cs 2 + o;|D ce cn ci ) 

( 8 ) 


Cutoff and Resonant Frequencies 

Explicit expressions for the resonant and cutoff frequencies may be readily 
obtained by setting k 2 equal to infinity and zero, respectively, in equa- 
tion (l). Cutoff frequencies are designated by subscripts beginning with o. 
Resonant frequencies are designated by the superscript oo. For each propagation 
angle there are usually three cutoffs and three resonances that are designated 
in the order of increasing frequency by the subscripts 1, 2, and 3. 

Cutoff . - For propagation at all angles except 0°, 

(“ci " ^ce) + V (“ci “ “ce) 2 + 4 Ke 0D ci + °>p) 

“ol = 2 * (9a ' 
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and 


“o3 


“Ki 



cd C 6 ) 2 + 4(cn ce o^ cl + co |) 
2 


(9b) 


%2 = ^ 


( 10 ) 


At 0°, co Q 2 does not exist. The numerical order of the subscripts in equations 
(9) and (10) prevails whenever C 0 p > cu ci or n > 0.03 B 2 . If co p < co ci , equa- 
tions (9) represent co 0 3 and oq-,2 and equation (10) represents cr> 0 ]_. 


Resonance . - At 6 = 0 °, 



(ID 


while at 9 = 90°, 
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(“ce + 
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J (12b) 


where and are referred to as the lower and upper hybrid frequencies, 

respectively. 


Resonant frequencies at an arbitrary propagation angle are determined by 
solving for the roots of the expression = 0 that was obtained by setting 
= oo in equation (l). Hence, 

3 2 r 

(1 + tan 2 ©) (co 2 ) - (l + tan 2 0)(co 2 e + cn 2 .^ + cu^)(cn 2 ) + (1 + tan 2 e)cD 2 .j_CD 2 e 

+ (^ ci cu ce )tan 2 0 + cd|(cu§ ± + co 2 e )J (cn 2 ) - a^a%±(D% e = 0 (13) 

This expression is also considered in reference 5. 


Phase and Group Velocities 

The phase and group velocities for the various modes are given by the ex- 
pressions 
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and 


U = 


CD 

k 


V = 

dk 

These are derived from equations (2), (3), (7), and (8). Only the expressions 
for group velocity are presented here. For the principal modes they are 

v = jiill ( 1 + — A . XI + JL <£ f 1 + 1 V i " 1 

1 [ 2 [_CD \0O - CD c i to + 03 ce yj 2k 7 c 2 yco - CD ci CD + CD ce Jj 


m Jirlz _ ( 1 

I 2 [CD *" \CD + 0 
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CD 


Ki 
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o|)cD 2 + (4 ± (4 e + o|cD ce CD ci 
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(14) 


For the principal modes, the group velocity is in the direction of wave propaga- 
tion. For waves propagating at some arbitrary angle to the magnetic-field di- 
rection, however, the group velocity is generally not in the direction of wave 
propagation. 

The expressions presented in the sections Dispersion Relation and Cutoff and 
Resonant Frequencies are sufficient to determine the major characteristics of 
electromagnetic waves propagating through a cold, collisionless plasma. There 
are several ways in which these equations may he simplified depending upon the 
approximations that are made. In appendix B three approximate treatments are 
discussed and compared to the Astrom model. These are (l) the microwave approxi- 
mation (ref. 6), (2) the Stix approximation (ref. 7), and (3) the simplified 
Ohm’s law (ref. 8). 


RESULTS AND DISCUSSIONS 

In the section THEORY, Astrom’ s dispersion relation was used to derive the 
characteristic wave equations in a plasma. These equations are valid for both 
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plane vaves and cylindrical waves propagating through a cold, collisionless 
plasma immersed in a uniform magnetic field. Except for figures 13 to 15, the 
data for all curves presented in this report were calculated from Astrom's equa- 
tions. Figures 13 to 15 were obtained from equations presented in appendix B. 
Calculations were performed for a fully ionized atomic hydrogen gas. Table I 
provides a summary of all the data covered by the curves of figures 1 to 18. 


Wave Number of Principal Modes 

Figure 1 consists of curves of the wave number squared as a function of the 
excitation frequency for the principal modes (P, S, R, and L) for a range of 
densities from 10 to lO^ 5 ions per cubic centimeter and at magnetic fields of 
10^, 10^, 10 and 10 5 gauss. A wavelength scale is also included. The dashed 
sections of the curves that appear in the figures are extrapolations from the 
calculated solid sections. Cutoff and resonant frequencies are indicated by 
triangular symbols at the bottom and top of the figure, respectively. Figures 
1(a) and (d) present the data for transverse wave propagation (0 = 90°) where 
the electric field of the wave is alined with the directions of the magnetic 
field and perpendicular to it, respectively? figures 1(b) and (c) present the 
data for longitudinal propagation (0=0°) where the electric field is perpendic- 
ular to the direction of the magnetic field and is left and right circularly 
polarized, respectively. 

Plots of k 2 as a function of co, such as those in figure 1, may be used 
directly to determine the feasibility of transmission experiments in which phase 
shift is measured. It is interesting to note that it is possible at certain con- 
ditions to determine an ion density with the use of the S mode, while it is 
impossible with the use of the P mode. For example, the maximum density that 
can be determined by a P mode at a frequency of 10ll radians per second is 
about 3xl0-*- 2 ions per cubic centimeter (fig. 1(a)). If the plasma were immersed 
in a magnetic field of 10^ gauss, densities up to 9x10-^ ions per cubic centi- 
meter could be determined with the use of the S mode (fig. l(d-3)). Conse- 
quently, by rotating the waveguide 90° about its axis, it would be possible to 
increase the upper density limit of the diagnostic system. Increasing the field 
strengths extends the upper density limits of the S mode to larger values. 
Another example of the use of the curves is illustrated in figure l(d-3). At a 
frequency of 10® radians per second, the validity of an experimental determina- 
tion of plasma properties with the S mode for a 10-centimeter-diameter plasma 
would be questionable for densities less than 10^ ions per cubic centimeter, 
because the wavelength of this mode is the same magnitude as the plasma size. 


Wave Number of Oblique Waves 

Wave numbers that have been presented up to now are for the principal modes 
in which propagation normal to the magnetic-field direction is uncoupled from the 
parallel propagation. For oblique waves, however, the components of the propa- 
gation vectors k x and k z , which are respectively normal to and parallel to 
B, are related to each other and to k through the relations 
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k2 = 4 + k z 

k x = k sin 0 
k z = k cos 0 

These x- and z-components of the propagation vector are related to the transverse 
and parallel rave numbers that appear in the analysis of bound plasmas (ref* 8)* 
These components of the oblique raves (along -with components of wavelengths) are 
presented in figures 2 to 8 as functions of o> for 9 = 0.1°, 1°, 10°, 45°, 80°, 
89°, and 89*9° and for representative values of B and n. The magnetic-field 
strengths covered by these plots are 10 2 , 10 3 , and 10^ gauss, and the particle 
densities are 10^ and lO-^- 3 ions per cubic centimeter* 

When the propagation angle is varied, the dispersion relation can be altered 
markedly in some regions of the frequency spectrum and remain relatively un- 
changed in other regions. In many experiments the exact propagation angle is 
difficult to determine, such as when waves are generated by induction coils. It 
becomes important, therefore, to assess the manner in which the wave number 
changes with the angle. 

Sketches (a) to (d) are presented to assist in the interpretation and use 
of figures 2 to 8 and to illustrate the effects of varying the propagation angle. 
In these sketches, the wave numbers of the principal modes are compared with 
components of the waves propagating at angles near 0° and 90°. 

Sketch (a) consists of curves taken from figures 1(b), 1(c), 4(a-2), 4(b-2), 
and 9(b). This sketch presents k 2 as a function of cd for the R and L 
modes and k z for a 10° propagation angle. As might be expected, k| at 10° is 

very similar to the wave numbers of the R and L modes over most of the fre- 
quency spectrum* There are, however, portions of the 10° dispersion relation, 
marked by crosshatching, that do not have counterparts In the R and L modes* 
This modification results in a resonance, which occurs at the plasma frequency, 
that cannot be predicted by the R and L modes. Near the 90° propagation 
angle the wave number might be expected to approximate the transverse S or P 
modes rather than the R and L modes* 

In sketch (b), which is taken from the curves of figures 1(b), 1(c), 6(a-2), 
6(b-2), and 9(b), k 2 for 80° propagation is compared to the R and L modes. 

Although the general shapes of the 80° curves bear some similarity to the R and 
L curves, the magnitudes differ appreciably. Obviously, the R and L modes 
may not be used to approximate the k z relation for large angles* 
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Sketch (c) is composed of curves taken from figures 1(a), 1(d), 6(b-l), 
7(b-l), 8(a), and 9(a). In this sketch k^ for propagation angles of 80°, 89°, 

and 89.9° is compared with the S and P wave numbers. Two interesting details 
should be noted in sketch (c). First, a resonance exists between (& c ± and zero 
frequency for propagation less than 90°; this resonance disappears at 90°. 


S mode 



8 , 

deg 

90 

89. 9 
89 
80 


(c) 


Second, for frequencies less than o^, the wave numbers change very rapidly as 
the angle is varied near 90°. It follows from sketch (c) that, for G less 
than 89.9°, k x cannot be approximated by the S and P wave numbers. In con- 
trast, in sketch (a) k| can be approximated by the R and L modes for angu- 
lar deviations as large as 10°. 

Another matter of passing interest is illustrated in sketch (c). For nega- 
tive values of k^, a branch of the P curve intersects a branch of the S 
curve. Their counterparts for angles less than 90° never intersect, however. 

As illustrated in sketch (d), there appears to be an interchange of branches as 
G approaches 90° (see branch marked (x)). 
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2 

x 


0 



(a) 

The trends for imaginary vave numbers are indicated in these sketches. The 
trends are important since the damping distance associated with wave phenomena 
can be represented by d = where ki m signifies the imaginary compo- 

nent of the wave number. In the absence of collisions, k is either purely real 
or purely imaginary when co is assumed real as in this report. To give an idea 
of the magnitude of the numbers involved, the specific case n - 10^-3 ions per 
cubic centimeter and B = 10 4 gauss is plotted in figure 9. The magnitude of 
ki m over most of the various branches readily explains why it is necessary to 
consider the imaginary region of the dispersion relation in wave propagation 
problems. As an example, the damping distance for k^ when Q = 90°, n = 10 13 

ions per cubic centimeter, and B = 10 4 gauss is roughly 0.17 centimeter (fig. 
9(a)). When 6 = 10°, however, the average damping distance is of the order of 
10 centimeters. If the characteristic length of an experimental plasma were 
much less than this latter damping distance, the plasma would be transparent over 
a large part of the imaginary- wave -number region. The plasma would then trans- 
mit a disturbance at both higher and lower frequencies than the transmission 
cutoff frequency, and no sharp and well defined cutoffs would exist. 
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In summarizing the discussions of sketches (a) to (d), it is worth noting 
the following for positive values of k 2 : Over the range of angles from 0° to 

80° the dispersion relations hear some degree of similarity to each other, while 
over the range 80° to 90° the low-frequency regions of the dispersion relation 
change in a drastic manner when the angle is slightly varied. This sensitivity 
to wave propagation direction near 90° is pointed out for the case of resonance 
in reference 5. 


Cutoff and Resonant Frequencies 

Cutoff and resonant frequencies are particularly interesting because of 
their relation to density measurement and power absorption in a plasma. At cut- 
off k 2 goes through zero and, therefore, is a dividing point between a propa- 
gation (k 2 > 0) and a nonpropagation (k 2 < 0) region. Similarly, at resonance 
k 2 goes from +oo to -oo and thereby again defines a separation between propa- 
gating and nonpropagating regions. In a real plasma at cutoff most of the signal 
is reflected by the plasma, while at resonance most of the signal will either be 
reflected or absorbed. In either case the signal can be attenuated appreciably 
as it passes through the plasma. If either the cutoff or resonant frequencies 
have a detectable variation with density, they may conceivably be applied to 
measure plasma density. 

Cutoff frequency is plotted against density in figure 16. The frequency 
<%p of figure 16(a) is the plasma frequency and is independent of magnetic 
field. The upper and lower hybrid resonant frequencies of the S mode are 
plotted in figure 17 as a function of ion density with the magnetic field as a 
parameter. The resonances of the R and L modes, which are the electron and 

ion-cyclotron frequencies, respectively, are plotted in figure 18 as a function 

of magn.etic-field strength. 

The resonant frequencies for oblique waves are plotted in figure 10 as a 

function of the propagation angle for magnetic-field strengths of 10 2 , 10^, and 

10^ gauss and for densities of 10-*-*- and 10^ ions per cubic centimeter. These 
frequencies are the roots of equation (13). Several interesting phenomena are 
illustrated in these curves. It is observed that the first resonant frequency 
of is practically constant between 0° and 80°. Above 80° the resonant frequency 

drops off very rapidly and goes to zero as 6 approaches 90°. This is consist- 
ent with the rapid change in the dispersion relation near 90° depicted in 
sketch (c) for the lowest resonance. The second resonance ajg continually de- 
creases as 9 goes from 0° to 90°. Therefore, in an experiment, if the propa- 
gation angle departs from that expected, the resonant frequency will shift from 
the expected value according to the presented curves. The third resonant fre- 
quency dj is fairly constant with angle, except, perhaps, at high field 
strengths (B > 10^ gauss) and large densities (n > 10-*-^ ions/cm^), where it in- 
creases with increasing angle. 

Another feature of the dispersion equation implied by the sketches and de- 
picted explicity in figure 10 is that the second resonant frequency approaches 
the smaller of either o>p or (D ce as 9 approaches 0°. The third resonance 
approaches the larger of the two frequencies. 
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Phase and Group Velocities of Principal Modes 


Phase velocity U is plotted as a function of cd in figure 11 for densi- 
ties of lO 1 ^ 10 13 , and 10)15 ions per cubic centimeter and for magnetic fields 
of 10 2 , 10 3 , 10^, and 10 5 gauss. A general property evident from all the plots 
is that the phase velocity for frequencies below the nonpropagating hand (or 
stop hand) is less than the speed of light, -while above the highest cutoff fre- 
quency it exceeds the speed of light. 


To complete the presentation, the group velocity V is plotted as a func- 
tion of co in figure 12 for densities of 10H, lO^- 3 , and 1015 ions per cubic 
centimeter and for magnetic fields of 10 2 , 10 3 , 10^, and 10 5 gauss. It may he 
noticed from the plots that the magnitude of the group velocity in the frequency 
region above the highest cutoff frequency is very near the speed of light and, 
consequently, would he difficult to measure experimentally within the presented 
frequency range. In the frequency region below the lowest transmission cutoff, 
however, the group velocity is small enough to he determined experimentally. 
Although the group velocities in the lower frequency regions of experimental in- 
terest are large, plasma densities could he determined by using waves of R or 
L mode if the propagation path length were sufficiently long to obtain a meas- 
urable time delay of the signal. As an example, for the plasma conditions pre- 
sented in figure 12(b-l) and a propagation length of 100 centimeters, a signal 
at a frequency of 10 7 radians per second would propagate through the plasma as 
an R mode with a time delay of the order of 2.5 or 25 microseconds, depending 
upon whether the density was 10^- 3 or 10^ 3 Ions per cubic centimeter, respec- 
tively. 


CONCLUDING REMARKS 

The characteristics of electromagnetic wave propagation through a cold, col- 
lisionless plasma immersed in a magnetic field have been presented in dimensional 
form so the experimentalist may have quantities in laboratory units. The wave 
numbers, wavelengths, phase velocity, and group velocity of the principal modes 
determined from Astrom ! s equation are presented as functions of frequency. Wave 
numbers (wavelengths) for oblique waves at various propagation angles are also 
presented as functions of frequency. These values are of interest in bound 
plasma experiments. In addition, the cutoff and resonant frequencies for both 
the principal and oblique waves are presented for a range of plasma conditions. 

In examining the results for oblique waves, it was noted that, over the range of 
propagation angles from about 80° to 90° with respect to the magnetic-field di- 
rection, the wave characteristics showed a marked sensitivity to propagation di- 
rection. This effect must be considered in the design of experiments and In the 
interpretation of experimental findings. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 7, 1963 
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S*YMBOLS 

dispersion relation coefficients 
approximate dispersion relation coefficient 

magnetic-field strength, gauss 
magnetic-field vector, gauss 
velocity of light, cm/sec 

damping distance 
electric field vector, esu 
electron charge, esu 
current vector* ahamp 

modified dielectric coefficient* Astrom ! s dielectric coefficient 
multiplied hy oy^/c^* K g = (K^ + ¥L^)/2 

wave number, 2rt/A* cm“^* 

propagation vector 

mass* g 

particle density, ions / cm^ 
phase velocity* oo/k* cm/sec 
group velocity* doo/dh* cm/sec 
mass velocity, cm/sec 
charge number of ions 

angle between direction of wave propagation and magnetic- field di- 
rection* deg 

wave length* cm 

wave frequency, radians/sec 

electron cyclotron frequency* Be/m e c* radians/sec 
ion cyclotron frequency* ZeB/m^c* radians/sec 



of 0 

OO 

OO 

^r 

“oZ 


“os 

“osl 

%s2 

“ol 

^02 

“o3 

% 

“pe 

“►pi 


resonance frequency of oblique waves 


resonance frequency 

of 

L mode, CD c j_, radians/sec 


resonance frequency 

of 

R mode, oo ce , radians/sec 


cutoff frequency of 

L 

mode, 

%i (eq- (9a)) or a> o2 , 

radians/sec 

cutoff frequency of 

P 

mode, 

%2 ( e< 3* ( 10 )) or 

radians/sec 

cutoff frequency of 

R 

mode, 

^03 ( e< l- ( 9*t> ) ) , radians/sec 

cutoff frequency of 

S 

mode, 

radians/sec 



lower cutoff frequency of S mode, eq. (9a), radians/sec 
upper cutoff frequency of S mode, eq. (9b), radians/sec 
lower cutoff frequency eq. (9a), radians/sec 
middle cutoff frequency, eq. (10), radians/sec 
upper cutoff frequency, eq. (9b), radians/sec 


plasma frequency, |4?tn e e 2 ^^- + m^)] ^ 9 radians/sec 
electron plasma frequency, (4rtn e e 2 /m e )"^ 2 , radians/s 


ion plasma frequency, I 4 j:Z 


00 

“si 

lower hybrid frequency of 

S 

mode, 

eq. (12a), radians/sec 

“s2 

upper hybrid frequency of 

s 

mode, 

eq. (12b), radians/sec 

“l 

approximate resonant frequency for 

S mode, eq. (B7), radians/sec 

“2 

approximate resonant frequency for 

S mode, eq. (B5), radians/sec 

C0 3 

CDci^ce + nadians/sec 




Subscripts: 




e 

electron 




i 

ion 
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I L mode, left circular polarization 

p P mode, axial polarization 

r R mode, right circular polarization 

s S mode, transverse polarization 

x rectangular coordinate, perpendicular to B 

z rectangular coordinate, parallel to B 
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APPENDIX B 


APPROXIMATE DISPERSION EQUATIONS 

The various approximations that are presented herein first appeared In the 
literature as developments of individual physical models. Because of their com- 
mon usage, wave numbers as functions of frequency were computed for these approx- 
imate dispersion relations and were compared with the results of Astrom f s expres- 
sions. The results of these comparisons are presented in this appendix. 


Microwave Approximation (ref. 6) 


When the ions are present only to maintain overall charge neutrality, a con- 
dition exists that is referred to as the microwave approximation. In this case 
the ions do not contribute to the total current, and an electron gas with charge 
neutrality results. Mower (ref. 6) has shown that the dielectric coefficients 
can be represented by the following equations (which can also be obtained from 
Astrorals expressions if it is assumed that » oo c j_ and If the inequality 
1 % « mj_ is used): 


k 


2 

l 


4 



4 



K 


s 



CO 2 - co§ e / 


Prom equations (Bl) and (B3), 


(Bla) 


(Bit) 


(B2) 


(B3) 



OXD, 


ce 


- o^e)(^ 2 + ^ce - ^pe) 
o?{o? - <x>§) 


(B4) 


where 


2 2 
“2 = 


ce + “pe 


(B5) 


As long as the frequency satisfies the inequality od » co c: £, the expression 
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developed for the microwave approximation yields practically the same results as 
Astrom’s expression; however, terms are dropped in this approximation that in- 
troduce errors in wave number whenever calculations are made near the cutoff 
(cD 0 p, co 0 j, co 0 g, and cd os2 ) or resonant (o^ 2 ) frequencies. These errors only be- 
come large when the transmission frequency is within a fraction of a percent of 
the cutoff or resonant conditions. 


Stix Approximation (ref. 7) 

Another model that is frequently used involves the simultaneous solution of 
Ohm’s law and the equation of motion of a fluid element. If the electron-ion 
mass ratio is not neglected relative to unity in the generalized Ohm's law for a 
collisionless plasma, results for Kj, Kp, K p , K g , and k 2 are obtained that 

are identical to Astrom’s results. In the calculations of reference 7 certain 
terms that contain the electron-ion mass ratio were neglected, and the wave num- 
bers are practically identical to Astrom’s results. Stix, in addition, used ap- 
proximate factors for a quadratic expression that appears in the equation for 
the dielectric coefficient. These approximations result in the expression 


where 


k 2 - — 
K s ~ o 

c c 


,2 (co 2 - oxocg 


- CO§) (l 


or + ^ce - “§) 


(co 2 - CO^) (go 2 - GO 2 ) 


v = “cki + Ki + 


-i 


pi 


J Ci/ 


(B6) 


(B7) 


<4 =4e + <4e ( B 5) 

In performing this reduction, it was necessary to neglect o ce 0D c i as com- 
pared to co^ • This approximation is equivalent to that which neglects the 
electron-ion mass ratio compared to unity. 

In a similar fashion it can "be shown that 



k 


2 

l 



(co - GD ci )(co + co ce ) 


4* 

(GO + CD cl )(C0 - CD ce ) 


(B8) 


(B9a) 


(B9b) 
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Since the only approximation made in the development of these equations was the 
neglecting of (& c ± compared to co ce , the error is large only in the vicinity of 
the frequencies that are within a fraction of a percent of the lower and upper 
hybrid resonances and the cutoff frequencies. 


Simplified Ohm l s Law (ref. 8) 


The previous approximation may be further simplified by neglecting the 
electron inertia term in Ohm T s law, that is, by using the equation 


E + X2<_5 

c 



X B 


0 


(BIO) 


in place of the generalized Ohm ! s law (eq. 2-12, ref. 1). This approximation is 
essentially equivalent to the condition that co « o> . For this case the diel- 
ectric coefficients are 
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(B12) 


The wave number for the S mode reduces to 


CO 2 

[^ce + (^ce^ci + cOp e )J [^ ce ” (^ce^ci + ^e) 

1 

C 2 

" “ci) 

" “pe“ci“ce^ 



For oblique waves the general dispersion relation (eq. (6)) becomes 


(B13) 


2 - K r )(k 2 - K 
K s k 2 - K r K z 

for this approximation. Note that Kp does not appear in the tan 2 0 relation. 

The component wave numbers may be written in the same form as equation (5), 
namely. 


(B14) 


tan 2 0 = 


* 
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sin 2 © 


k !l = 


* 1 + 4 * 1 - 4 -^i y i 


2sd-\ 


k x2 ~ 


*i - V s ! - iJ *i Y i . 


sin 2 © 


k zl 


__ % + V^i ~ ^ 1 2 


2 j/i 


cos © 


2 _ a i ~ V^i ~ 4j ^i 


z2 


2 s/-, 


^ cos 2 © 


(B15) 


The expressions for and *i have been modified, however, so that in 

this case 


s/-^ = cos 2 © 


= K s (l + cos 2 ©) > 


*L = K rK 2 


(B16) 


These expressions are identical to the dispersion -relation presented in refer- 
ence 8 if the proper substitutions for k and Q are made* 


Figure 13 presents curves calculated from the dispersion relation obtained 
by this method. When particle currents are negligible compared to displacement 
currents. Ohm's lav does not affect the results, and, consequently, the restric- 
tion cd « co ce is not required. Therefore, there are regions of the frequency 

spectrum, other than cd « CD ce , vhere the approximate and exact expressions 
agree. The plots of and k| for large magnetic fields (B = 10 5 , figs. 

13(T>-2), 13(c-2), l( c-4), and l(d-4)) demonstrate that agreement may exist for 
nearly the complete frequency range. Discrepancies betveen the tvo models for 
the example occur in narrov frequency bands centered around the upper resonance 
points (near <r> ce ), which, incidentally, are not predicted by the approximate 
relations. As the magnetic field decreases, there is a wider range of frequen- 
cies for which the approximate equation is unsatisfactory. The regions of inac- 
curacy are particularly large for low densities. These latter conclusions are 
exemplified by comparing figures l(b-l), l(c-l), and l(d-l) with 13(a-l), 

13(b-l), and 13(c-l), respectively. 

Figures 14 and 15 present the approximate dispersion relation for oblique 
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waves propagating at angles of 10° and 80°, respectively, as calculated for the 
simplified 0hm f s law approximation. Curves presenting the corresponding Astrom 
data are also included in the figures to illustrate, hy comparison, in what fre- 
quency range the results of this approximation can he applied. 
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TABLE I. - SUMMARY OF DATA PRESENTED IN FIGURES 1 TO 18 


Figure 

Curve description 

Model 

Wave propaga- 
tion angle, 

e, 

deg 

Magnetic- 

field 

strength, 

B, 

gauss 

Particle 

density, 

n, 

ions/cm 3 

1 

P 

k^ against cn 
(k 2 > 0) 

Astrom 

0, 90 

10 2 to 10 5 

10 10 to 10 15 

2 to 8 



2 2 

k x , k z against to 
(k 2 >0) 

Astrom 

0. 1, 1, 10, 45, 
80,89,89. 9 

10 2 to 10 4 

10 11 , 10 13 

9 

k 2 , k 2 against cd 
( k 2 < 0) 

Astrom 

1, 10, 45, 80, 
89,89. 9, 90 

10 4 

10 13 

10 

of 0 against 6 

Astrom 

0 to 90 

10 2 to 10 4 

10 11 , 10 13 

11 

U against cd 
(k 2 >0) 

Astrom 

0, 90 

10 2 to 10 5 

10 11 to 10 15 

12 

V against oo 
(k 2 >0) 

Astrom 

0, 90 

10 2 to 10 5 

10 11 to 10 15 

13 

k 2 against cd 
(k 2 >0) 

Simpli- 
fied 
Ohm 1 s 
law 

0, 90 

10 2 , 10 5 

10 10 to 10 15 

“i 

H 

V* 

H 

CJ1 

k 2 , k| against cd. 
(k 2 > 0) 

Simpli- 
fied 
Ohm 1 s 
law 

10, 80 

10 2 , 10 4 

10 11 , 10 13 

16 to 18 

cd 0 ,cd°° against n 

Astrom 

0, 90 

10 2 to 10 5 

10 10 to 10 15 
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Figure 1. - Wave number as function of frequency. Wave propagation angles, 0° and 90°; Astrom's model. 

Figure 2. - Wave number as function of frequency. Wave propagation angle, 0.1°; magnetic-field strength, 10 4 gauss; particle den- 
sity, lu 13 ions per cubic centimeter; Astrom's model. 

Figure 3. - Wave number as function of frequency. Wave propagation angle, 1°; Astrom's model. 

Figure 4. - Wave number as function of frequency. Wave propagation angle, 10°; Astrom's model. 

Figure 3. - Wave number as function of frequency. Wave propagation angle, 45°; Astrom's model. 

Figure 6. - Wave number as function of frequency. Wave propagation angle, 80°; Astrom's model. 

Figure 7. - Wave number as function of frequency. Wave propagation angle, 89°; Astrom’s model. 

Figure 8. - Wave number as function of frequency. Wave propagation angle, 89.9°; magnetic-field strength, I0 4 gauss; particle den- 
sity, 10 13 ions per cubic centimeter; Astrom's model. 

Figure 9. - Wave numbers in imaginary {nonpropagating) region as functions of frequency for various oblique angles. Magnetic-field 
strength, 10 4 gauss; particle density, 10 13 ions per cubic centimeter. 

Figure 10. - Resonant frequency as function of oblique wave propagation angle. Astrom's model. 

Figure 11. - Phase velocity as function of wave frequency. Astrom's model. 

Figure 12. - Group velocity as function of frequency. Astrom's model. 

Figure 13. - Wave number as function of frequency. Wave propagation angles, 0° and 90°; simplified Ohm's law approximation. 

Figure 14. - Wave number as function of frequency. Wave propagation angle, 10°; simplified Ohm's law. 

Figure 13. - Wave number as function of frequency. Wave propagation angle, 80°; simplified Ohm's law. 

Figure 16. - Cutoff frequencies as functions of density. 

Figure 17. - Resonance frequencies of principal modes as functions of particle density. 


Figure 18. - Resonance frequencies of R and L modes as functions of magnetic-field strength. 




Wave frequency, 
(a) P mode; wave 
Figure 1. - Wave number as function of frequency. 
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to, radians/a ec 
propagation angle, 90°. 

Wave propagation angles, 0° and 90°; Astrom’s model. 
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Figure 1. - Continued. Wave number as function of 
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Wavelength. 



(d) Continued. S mode; 
Figure 1. - Continued. Wave number as function of 
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Figure 1. 


(d-4) Magnetic-field 
(d) Concluded. S mode; 
- Concluded. Wave number as function of 




Wave frequency, o>, radians/sec 
(a) x-component. 


Figure 2. - Wave number as function of frequency. Wave propagation angle, 0.1°; magnetic-field strengths 
10 gauss; particle density, 10- ^ ions per cubic centimeter; Astrom's model. 
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Figure 3. - Wave number as function of frequency. Wave propagation angle, 1°: Astrom's model 
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Figure 3. - Continued. Wave number as function of frequency. Wave propagation angle, 1°: Astrom’s model 
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Figure 3. - Concluded. Wave number as function of frequency. Wave propagation angle, 1°; 
Astrom’s model. 
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Figure 4. - Wave number as function of frequency. Wave propagation angle, 10°; Astrom's model. 
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Figure 4. - Continued. Wave number as function of frequency. Wave propagation angle, 10°; 
CD Astrom r s model. 
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Figure 4. - Continued. Wave number as function of frequency. Wave propagation angle, 10°: Astrom's model. 
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Figure 4. - Concluded. Wave number as function of frequency. Wave propagation angle, 10°; 
Astrom’s model. 
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Figure 5. - Concluded. Wave number as function of frequency. Wave propagation angle, 45°; Astrom’s model. 
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Figure 6. - Wave number as function of frequency. Wave propagat'on angle 
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Figure 6. - Continued. Wave number as function of frequency. Wave propagation angle, 80°; Astrom's model 
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Figure 6. - Continued. Wave number as function of frequency. Wave propagation angle, 80°; Astrom's model. 
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Figure 7. - Continued. Wave number as function of frequency. Wave propagation angle, 89°; Astrom’s model. 
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Figure 8. - Wave number as function of frequency. Wave propagation angle, 89. 9° j magnetic-field strength, 10 4 gauss; particle density 
10 13 ions per cubic centimeter: Astrom's model. 
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Figure 8. - Concluded. Wave number aa function of frequency. Wave propagation angle, 89.9°; magnetic-field strength, 
10 gauss; particle density, 10^3 ions per cubic centimeter; Astrom's model. 






Wave propagation angle, 0, deg 

(b) Magnetic-field strength, 10^ gauss. 

Figure 10. - Continued. Resonant frequency as function of ob 
lique wave propagation angle. Astrom's model. 
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Figure 11. - Continued. Phase velocity as function of wave frequency. Astrom's model 
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Figure 11. - Continued. Phase velocity as function of wave frequency. Astrom's model 
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Figure 11. - Concluded. Phase velocity as function of wave frequency. Astrom’s model 
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Figure 12. - Continued. Group velocity as function of frequency. Astrom's model 
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Figure 12. - Concluded. Group velocity as function of frequency. Astrom's model 
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Figure 13. - Wave number as function of frequency. Wave propagation angles, 0° and 90°; simplified Ohm’s law approximation. 
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Figure 13. - Continued. Wave number as function of frequency. Wave propagation 
angles, 0° and 90°; simplified Ohm's law approximation. 
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Figure 13 


Continued. Wave number as function of frequency. Wave propagation angles, 0° and 90°; simplified Ohm's law approximation 
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Figure 13. - Continued. Wave number as function of frequency. Wave propagation angles, 0° and 90°; simplified Ohm's law approximation 




Wavelength 


H 

O 

CD 



Wave frequency, o>, radians/sec 
(c-l) Magnetic-field strength, 10 2 gauss, 
(c) S mode; wave propagation angle, 90°. 


Figure 13. - Continued. Wave number as function of frequency. Wave propagation angles, 0° and 
90 5 simplified Ohm’s law approximation. 
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Figure 15. - Continued. Wave number as function of frequency. Wave propagation angle, 80°; simplified Ohm's law. 
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(b) Magnetic -field strength, 10^ gauss; particle density, 10^3 ions per cubic centimeter. 

Continued. Wave number as function of frequency. Wave propagation angle, 80°; simplified 0hm l s law. 
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Figure 15. - Concluded. Wave number as function of frequency. Wave propagation angle, 80°; simplified Ohm^s law. 
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Figure 16, - Continued. Cutoff frequencies as functions of density 
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Figure 16. - Concluded. Cutoff frequencies as functions of density. 
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Figure 17. - Resonance frequencies of principal modes as functions of particle 
density. 
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Figure 17. - Concluded. Resonance frequencies of principal inodes as functions of 
particle density. 
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